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A new class of solutions to Laplacian growth (LG) with zero surface tension is presented and 
shown to contain all other known solutions as special or limiting cases. These solutions, which are 
time-dependent conformal maps with branch cuts inside the unit circle, are governed by a nonlinear 
integral equation and describe oil fjords with non-parallel walls in viscous fingering experiments in 
Hele-Shaw cells. Integrals of motion for the multi-cut LG solutions in terms of singularities of the 
Schwarz function are found, and the dynamics of densities (jumps) on the cuts are derived. The 
subclass of these solutions with linear Cauchy densities on the cuts of the Schwarz function is of 
particular interest, because in this case the integral equation for the conformal map becomes linear. 
These solutions can also be of physical importance by representing oil/air interfaces, which form oil 
fjords with a constant opening angle, in accordance with recent experiments in a Hele-shaw cell. 
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I. INTRODUCTION 

Background and motivation: The Laplacian growth 
(LG) is a free boundary motion governed by gradient 
of a harmonic field. It describes numerous physical pro- 
cesses [l[ with a moving boundary between two immis- 
cible phases (an interface) far from equilibrium and pro- 
foundly interconnects various branches of mathematical 
physics 0, [H • LG in two dimensions is of special impor- 
tance since in the zero surface tension limit it has a re- 
markably rich integrable structure [1,0, @] and an impres- 
sive list of exact solutions in the form of conformal maps 
with pole @, S, Q and logarithmic [H 03, El, El time- 
dependent singularities (many of these solutions cease 
to exist in a finite time due to the interface instabil- 
ity.) In addition, there exist a family of explicit solu- 
tions [3, EH, EH with time-independent branch cuts with 
fixed endpoints which describe growth of self-similar air 
'fingers' in a wedge geometry (or in the plane with the 
rotational Zjy-symmetry imposed). 

Do other solutions exist, apart from those mentioned 
above? Answering 'yes' to this question, we hereby 
present a new family of conformal maps with time- 
dependent cuts, study their dynamics and provide ex- 
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amplcs. This family is the most general solution to the 
LG problem (in the absence of surface tension) which 
can be written in terms of analytic functions. All ratio- 
nal, logarithmic and self-similar solutions known so far 
(and mentioned above) are shown to be special or limit- 
ing cases of these newly found ones. 

In this paper we consider a simply connected exterior 
2D LG problem with a source/sink at infinity. In the 
context of Hele-Shaw flows [l7| (the prototype of all 2D 
LG processes [13), it corresponds to a finite bubble of 
an inviscid fluid (air) surrounded by a viscous fluid (oil) 
with a source or sink at infinity. In what follows we hold 
to this hydrodynamic interpretation. 

The standard formulation of the exterior LG problem 
in 2D planar geometry is as follows: 



v = — Vp 
Ap = 
p = const 

V n = -0 n p 

. P(x,y) = -^rlogi?, when R = y/ ' x 2 + y 2 -> oo. 



in D(t) 

at r(t) 



(1) 



Here t is time, T(t) is the boundary (an oil/air interface) 
of an infinite planar domain D(t) containing infinity and 
filled with oil, p is the pressure field, v is velocity of oil, 
V n is normal velocity of the oil/air interface, d n stands for 
normal derivative at the interface, and Q is a pumping 
rate assumed to be time-independent (positive for the 
suction problem and negative for the injection one). 

The first equation in the system ([I]) is the Darcy law 
in scaled units, which determines the local velocity of oil 
in the Hele-Shaw cell. The second equation follows from 
continuity and incomprcssibility of oil, divv = 0. The 
third equation results from the fact that in fluids with 
negligible viscosity pressure is the same everywhere in 
the fluid, p = const, so, neglecting surface tension, one 
can set p to be the same constant along both sides of the 
interface. The fourth equation states that normal veloc- 
ity of fluid at the interface and that of the interface itself 
coincide (continuity). The last equation gives the asymp- 
totic of p far away from the interface in the presence of 
a single sink of strength Q placed at infinity. 

A reformulation of the LG problem in terms of complex 
variables and analytic functions is particularly instruc- 
tive. We conclude the introduction by a brief review of 
the time-dependent conformal map approach. 

The conformal map formulation. Let us introduce a 
time-dependent conformal map, z = f(w,t), from the 
exterior of the unit circle, \w\ > 1, in an auxiliary math- 
ematical w -plane to the domain D(t) occupied by oil in 
the physical plane z = x + iy. Then the interface at 
time t is a closed curve swept by f(e l< ^, t) as <j> runs from 
to 2tt. Below we sometimes do not indicate the time 
dependence explicitly and write simply f(w). It is con- 
venient to normalize the conformal map by the condition 
that oo is mapped to oo and the derivative at infinity is 
a real positive number r (the conformal radius), so that 



the Laurels expansion at infinity has the form 

f(w) = rw + u + ui/w + ui/w + . . . (2) 

(the coefficients v,i are in general complex numbers). 

It is well known [l9|, [2(| that the LG dynamics of 
the interface is equivalent to the following equation for 
/(e* t): 



(3) 



where subscripts stand for partial derivatives and a bar 
for complex conjugation. This nonlinear equation is re- 
markable because it possesses an infinite number of con- 
servation laws Q and an impressive list of non-trivial so- 
lutions with moving singularities 0, H, 0, E3, EH, E2, EH , 
which arc cither poles or logarithmic branch points ly- 
ing strictly inside the unit circle. In all these cases the 
conformal map f(w) thus admits an analytic continua- 
tion across the unit circle, so the function f(w) is actually 
analytic not only in its exterior but in some larger infinite 
domain containing it. Assuming that f(w,t) is analyt- 
ically extendable across the unit circle for all t in some 
time interval, one can analytical ly continu e the LG equa- 
tion ^ itself. Set f(l/w,t) = f(l/w,t), then equation 
([3]) can be rewritten as 

d w f(w, t)dj(l/w, t) - dtf(w, t)d w f(l/w, t) = 1/w (4) 

(we have set Q — tt that means a resealing of time units). 
Under our assumption the functions f(w, t) and /(1/w, t) 
have a common domain of analyticity containing the unit 
circle \w\ = 1, and w in equation (0| is supposed to be- 
long to this domain. In fact, for rational and logarithmic 
solutions, this equation holds everywhere in the w-plane 
except for poles and branch cuts of f(w,t) and f(l/w,t). 
We note in advance that the general multi-cut solutions 
constructed below have the same property. Namely, the 
functions f(w,t) and f(l/w,t) have a common domain 
of analyticity where they actually solve the analytically 
continued LG equation (U]) and not just ([3]). 

As said earlier, there are also self-similar Zjv- 
symmetric solutions 0, E3, EH, EH containing a time- 
independent fractional power singularity which corre- 
sponds to the interface self-intersection point at the ori- 
gin. In this case the analytic continuation is possible 
through every point on the unit circle except those that 
are mapped to the origin. Again, the functions f(w, t) 
and f(l/w,t) have a non-empty common domain of an- 
alyticity and equation ([4J is valid everywhere in this do- 
main. 



II. THE SCHWARZ FUNCTION AND 
INTEGRALS OF MOTION 

It is known @, H i, EI E3, HI that the partial dif- 
ferential equation ^ can be integrated in terms of the 
Schwarz function. It is an analytic function, S, which 
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connects f(w) and f(l/w) in their common domain of an- 
alyticity: /(l/u>) = S(f(w)). Equivalently, the Schwarz 
function for a curve T is an analytic function S(z) such 
that z = S(z) for z G T [13]. In other words, S(z) is 
an analytic continuation of the function z away from the 
curve. For analytic curves this function is known to be 
well defined in a strip-like neighborhood of the curve. If 
the curve depends on time, so does its Schwarz function, 
S = S(z,t). 

In terms of the Schwarz function the (analytically con- 
tinued) LG equation (g]) reads [21 1 



S t = d z log w 



(5) 



where w — w(z,t) is the function inverse to the f(w,t). 
Since d z logw(z) is analytic everywhere in the oil domain 
D(t) 1 both sides of equation © are free of singularities 
there. Therefore, all singularities of S(z,t) located in 
D{t) are time-independent. Equivalently, the function 



S+{z) 



s(0 



2ni 



(0) 



defined for z outside D(t) by the integral of Cauchy type 
(where the integration contour T has the standard anti- 
clockwise orientation) and analytically continued to D(t) 
is constant in time: d t S + (z,t) = 0. This implies conser- 
vation of harmonic moments tk Q which are coefficients 
of the Taylor expansion of S+(z) at z = 0: 



t k = 
1 

irk 



1 



2nik 



r(t) 
dx dy 



zdz 



k = 1,2,3,. 



D(t) 



(7) 



where we assume, without loss of generality, that the con- 
tour T(t) encircles the origin. In particular, t± = S+(0). 
Besides, it is easy to see that the physical time t is equal 
to the area of the air bubble divided by 7r: 



t 



r(t) 



z dz 



(8) 



Thus the LG dynamics is a process such that the area sur- 
rounded by the interface grows linearly with time while 
the harmonic moments ([7]) (or the S^-part of the Schwarz 
function ©) are kept constant. It is important to note 
that the function 5+ (i.e., the infinite set of harmonic 
moments tk) together with the variable t fix the whole 
Schwarz function S uniquely, at least locally in the vari- 
ety of closed analytic contours in the plane. In a certain 
sense (made more precise in [Hj]), these data serve as 
local coordinates in the infinite dimensional variety of 
contours. 

It is also worthwhile to mention that the function 
SL(z) = S(z) — S+(z) is analytic in D(t) and vanishes at 
infinity as 0(l/z) with the residue 



[S-{z)dz\ = —t 



(9) 



(the residue at infinity is defined as res^dz / z\ = — 1). 
The decomposition of the Schwarz function S = S+ + 5_ 
appears to be very useful. As it follows from properties 
of Cauchy-type integrals, the function S-(z) is given by 
the same integral ((6]) (with the opposite sign), where z 
now belongs to D(t). 

The conformal map and the Schwarz function are con- 
nected in the following way: 



S(z) = f(l/w) 
z = f(w) 



(10) 



These formulas make it clear that if the functions f(w) 
and f(l/w) have a common domain of analyticity con- 
taining the unit circle, then the Schwarz function is well- 
defined. They also imply the one-to-one correspondence 
between singularities of the function f(w) inside the unit 
disk and singularities of the function S(z) in D(t) (which, 
by the properties of the Cauchy type integrals, are the 
same as singularities of the function S+(z)). 

Let us illustrate the one-to-one correspondence be- 
tween singularities of f(w) and S+(z) by the example of 
poles in which case it is especially transparent. Suppose 
that the conformal map f(w) has a pole of order k at a 
point a inside the unit disk, so that the local behavior of 
/ near a is 



A 



(w — a) h 



less singular terms. 



(11) 



Then the function S = f(l/w) has a pole of the same 
order at the point w = I /a: 



S = f(l/w) 



A 



(w 1 — a) k 



less singular terms. (12) 



Because the point 1/a together with its small neighbor- 
hood lies outside the unit disk (where / is conformal), 
we can linearize the denominator in (|12p near z = /(1/a) 
and obtain 



S(z) = A 



r - /(1/a) 



1 A' 



less singular terms (13) 



(here / s (l/a) = d w f(l/w) at w = a). Thus S(z) neces- 
sarily has a singularity at z = /(1/a) of the same type 
as f{w) has at w = a. 

Now we are ready to outline the general strategy of 
integration of the LG problem. Taken initial data, i.e., a 
conformal map f(w) at t = to, one finds 5+ from ^ by 
rewriting it as integral over the unit circle in the mathe- 
matical plane: 



S+(z) 



1 

2~7Ti 



f(l/w) df(w) 
f(w) - z 



(14) 



In order to obtain the LG dynamics, one then should 
solve the inverse problem: given the S+(z) and time t, 
to recover the conformal map f(w,t). The last step is a 
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part of the inverse potential problem and thus is hard to 
implement in general. Formally, this is equivalent to a 
nonlinear integral equation which is easy to derive using 
the definition of the Schwarz function in the form /(C) = 
Integrating both sides of this equality with 
the Cauchy kernel 1/ (w — £) over the unit circle assuming 
that \w\ > 1, and using the fact that the ^--part of the 
Schwarz function does not contribute to the integral, we 
get the integral equation 



f(w) = rw + u + 



1 

2tt7 



CI 



5+(/(i/C))dC 



\w\ > 1. (15) 



The same procedure at \w\ < 1 gives a "complimentary" 
equation 



rw + uq 



1 

2ni 



5+(/(l/C))dC 



IC 1=1 



£ — w 



+ S-(f(l/w)), \w\<l. 

(16) 

As their direct consequence, we can write the following 
useful formulae for the coefficients of the Laurels series 



1 

2Tri 



S+tfil/Q)?- 1 ^, j>0, (17) 



I CI =i 

1 

2Tri 



S + (/(l/0)C" 2 dC + - 

ICI=i r 



(18) 



Substituting them back into ^ and using ([15]) , we get a 
set of integral equations 



/O) 



2iri 



rw 



u i w J + 

3=0 

C d g + (/(VC)K 

w d (w — C) ' 



\w\ > 1. (19) 



for any d = —1,0,1,2,.... Being taken together with 
(|17p . all of them are equivalent to the original equation 
(|15p which is reproduced at d = 0. However, depend- 
ing on the type of singularities of the function S + , one 
or another form may be more convenient than others. 
The choice d = — 1 gives an equation where none of the 
coefficients Uj enter explicitly: 



1 

2ni 



ICI 



wg + (/(i/Q)rfC 



|tu| > 1. (20) 



This equation combined with relation (|18|) accomplishes 
integration of the LG problem. Let f(w) = f(w\r) be 
a solution to equation (|20[) depending on the parameter 
r, then the time dependence is found from (|18[) which 
determines r as an implicit function of t: 



t = r z 



(21) 



Alternatively, one may use the relation 



t = 



1 

2ni 



\w\ = l 



f(l/w)df(w) 



(22) 



leading to the same results. 

In principle, this scheme provides a general solution to 
the LG problem but in a rather implicit form. However, 
for several important classes of functions S+(z) it can be 
made more explicit. In all effectively solvable cases, the 
r.h.s. of the integral equation simplifies drastically after 
shrinking the integration contour to singularities of the 
function <S + (/(l/£)) inside the unit disk. Some examples 
arc given in the next section. 



III. RATIONAL AND LOGARITHMIC 
SOLUTIONS 

In this section we apply the general method based on 
equations (|20|) , (|21|) to construction of LG solutions with 
poles and logarithmic singularities. 



A. Rational Solutions 

By rational solutions we mean solutions to equation Q 
whose only singularities are poles inside the unit disk (not 
mentioning a simple pole at infinity). The very fact that 
the rational ansatc is consistent with the LG equation, 
i.e., that the number of poles is conserved and singulari- 
ties of other types are not generated, is a consequence of 
the integral equation (fTS"!) . 

Let us take S+ to be a rational function of a general 
form 



A' 



T, 



K a 



' ^ (z 

1 fc=l v 



71 



0,0 



fe=l 



A: Z 



(23) 



where all the poles z m are in D(t) and T m ^ are arbitrary 
complex constants. Then the function has 
a pole at of order Kq and poles of orders K m at the 
points a m such that z m = f(l/a m ) and does not have 
other singularities inside the unit disk. For convenience, 
we set gq = 0. Therefore, the integral in the r.h.s. of (TT5|) 
is equal to the sum of residues at these poles. Calculating 
the residues, we obtain an expression for the conformal 
map in the form of a rational function of w with poles of 
orders not higher than K m at the points a m inside the 
unit disk: 



N K„ 



A m , k (t) 



m— k= 



\ {w - a m {t)) k ' 



(24) 



The pole at ao = is distinguished among the others 
because it does not move. All other poles as well as 
all coefficients depend on time. The coefficients A m ^ 
are expressed through derivatives of the function f(l/w) 
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at w — a m . The general formulae are quite compli- 
cated. In fact what we really need are inverse formu- 
lae which express the integrals of motion T m ^ through 
time-dependent coefficients of f(w,t). Obviously, 



T 



rn , k 



[{z-z m ) k - l S+(z)dz\ , m^O. 



then the expressions for T m .i and t acquire an especially 
compact and explicit form: 

T m ,i = A mA (t)d w f(l/w,t)\ w=!imit) = 

V a mW , (1 - ai(t)a m (t)) 2 ) 



and the function S+(z) here can be substituted for the 
S(z) because S-(z) is regular in D(t). Passing then to 
the mathematical plane, we obtain the full list of con- 
stants of motion: 

z m = f(l/a m (t),t), m^O (25) 
T mi to,k = res s -i (t) [(f(w,t)-z m ) k - 1 f(w-]t)df(w,t)] 

T , k = -re Soo [U{w,t))- k - l f(l/w,t)df{w,t)\ . 

We note that T m ^ have the meaning of "Whitham times" 
for a general Whitham hierarchy which covers the LG 
problem with zero surface tension. Each Whitham time 
T m ,k is "coupled" to its own type of singularity. From 
this point of view, formulae ([23)) are relations of the hodo- 
graph type which provide a solution of the full genus zero 
Whitham hierarchy in an implicit form (24|. 

Algebraic equations (|25p express time-dependent pa- 
rameters r(t), uo(t), a m (t) and A m ^{t) implicitly via the 
constants of motion z m and T mk . To obtain a complete 
set of relations, we also need the equation (|22|) containing 
time t explicitly. In our case the integral is reduced to a 
sum of residues: 



-reSoo[/(w 1 ,t)rf/(w,^)]-^res a -i (t) [/(w 1 ,t)df(w,t)} 



t = 



Comparing with (|25[) at k = 1, we see that the residues 
at l/d m are just constants T m i. Therefore, 



t 



,[f(l/w,t)df(w,t)} 



(26) 



So, we have N + K + 1 complex algebraic equations, 
where K = ^2^ =0 K m , and one real algebraic equation 
for N+K+l complex parameters and one real parameter 
which determine them as functions of time and N+K+l 
constants of motion. Given initial conditions in the form 
(|2"3")) , a solution to this system allows one to recover the 
time-dependent conformal map, according to the strategy 
outlined at the end of Section [TTJ 

In some important particular cases the general expres- 
sions given above become simpler and more explicit. For 
example, if none of the poles a m lies at the origin (i.e., 
A ok = for all k > 1), then T k = for k > 1, 
To,o = /(0,i) and 



t 



-(t)f(o,t)-E T » 



(27) 



(here f'(0,t) is the derivative of f(w,t) at w = 0). If, 
moreover, all poles in (|24p arc simple and none of them 
lies at the origin (K m = 1 for m ^ and Kq = 0), 



m 

Ai,i{t)A m ,i{t) 



I.) 1 1 



(1 - ai(t)d m (t)) 2 ' 



(29) 



B. Logarithmic Solutions 

For the class of logarithmic solutions, the function S+ 
is taken in the form 



JY 



S + {z) = h + J2A m log 1 , 



(30) 



where all the branch points z m are in D{t) and t\, A m 
are arbitrary complex constants (note that this t\ is the 
first harmonic moment from {7}). If 



N 



A 



J2 A ™ 



then there is also a branching at infinity. This function is 
multi- valued and one should fix a single- valued branch. 
In a neighborhood of the origin, we define it by the con- 
dition that >!?-|_(0) = t\. In order to continue it unam- 
biguously to D(t), it is necessary to introduce a system 
of cuts connecting all the branch points in such a way 
that the domain D(t) \ {all cuts} be simply connected. 
There are many ways to draw the cuts. Although all of 
them are ultimately equivalent, it would be convenient 
for us to meet some requirements natural for the growth 
problem: to respect the democracy among the branch 
points z m and to ensure that the cuts always remain in 
D(t). For the former, let us fix a point q in D(t) and 
make cuts from q to all z rn (and, if necessary, to infin- 
ity). For the latter, it is natural to choose q = oo since 
oo is the only point in D(t) which remains there forever 
irrespectively of initial conditions. So, we fix a system of 
(non-intersecting) cuts from the points z m to oo. Note 
that this choice of cuts implies that S+(z) does not have 
a definite value at z = oo, even if oo is a regular point, 
because the limit depends on a particular way of tending 
z — > oo. 

Reconstructing / from (|15[) in a similar way as for ra- 
tional solutions, we introduce the points a m such that 
f(l/a m ) = z rn inside the unit disk and notice that the 
integral in (Tl?)) is shrunk to residues of the differential 

N 

dS+(f (1/0) = ]T A m dlog(f(l/0 - f(l/a m )) 
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which are easy to calculate. In this way one obtains 

JV 

f(w, t) = r(t)w + u„(t) + ^ A m \og(w - a m (t)), (31) 



JV 



A m = 0. 



m— 

Again, we distinguish a possible branch point at the ori- 
gin and set do = 0. By construction, the cuts of the 
function f(w) go from to a rn inside the unit disk. Sim- 
ilarly to the Schwarz function, f(w) with this choice of 
cuts does not have a definite value at w = 0, even if is a 
regular point, because the limit depends on a particular 
way of tending w — > 0. 

The constants of motion are given by 



iio 



1 

2~7ri 

JV 



f(l/w)df(w) 



+ A m log 



/O) 

(2 m Z. 



r 



m— 1 

2m = f(l/a m (t),t) 

X 

'/am(t) + u{t) +Y. Al - a ' (*)am (*)), (32) 



2=1 



in agreement with previous results [lCJ, [ll|, ll_2j, ll_3j . For- 
mulae (|2Tj) and ([22)) allow one to derive few differently 
looking equivalent expressions for t. The simplest one 
reads 



t = r 2 + ^ A r 



\z m -u - ra n 



(33) 



m=l 



The others can be also useful in some situations: 

t = rd w (/» + A Q log f(l/w))\ w=Q + A m z m , (34) 



t = r 2 + ^ A[A m log(l - aia m ). 



(35) 



The solution contains N + 1 complex (a m and Mo) and 
one real (r) time-dependent parameters. They are to be 
determined from the system of N + 1 complex equations 
(|32| and one real equation (|34| . 



C. "Mixed" Rational-Logarithmic Solutions 



where a single-valued branch is fixed by the condition 
that 



N K m 
m=lfc=l m 



0,0 



(37) 



The equation (fTS]) gives 



/(«, t) = r(t)w + Uo (t)+J2J2 ( AmA t^ 



N 



22 A m l0g(w - a m (t)) 



JV 

E 

m=0 



Am . 



(38) 



with the same convention ao = 0. The constants of mo- 
tion z m , Tm.fc arc expressed through the time-dependent 
parameters of the conformal map by means of formulas 
similar to (|25|) . In fact one can represent the "hodograph 
relations" (f2"S")) in a form which is suitable for logarithmic 
and mixed cases as well: 

z m = f(l/a m (t),t) 

T m ^o,k = - ^resi/ 5m ( t ) [(f(w,t) - z m ) k df(l/w, t)] , 

ZWo = - ^re Soo [(f(w, t))- k df(l/w, t)] , (39) 

The residue at infinity is well defined since the differential 
df(l/w,t) is single- valued. Note that the coefficient A m 
can be formally understood as T m fl\ A m = —kT m ^ at 
k = 0. An expression for To,o follows from (|3T[) and the 
integral formula for t\ in (f3"2"| valid in all cases: 



Tq,o — 



1 

2iri 



JV K„ 

I, ^gp -EEy- (40) 

H=l m=lfe=l m 

Different equivalent versions of the formula for t can be 
obtained from (J2TJ) or ([22]) . 

Finally, we note that rational terms in f(w) can be re- 
garded as a special (singular) limiting case of logarithmic 
ones with merging branch points. For instance, 

i- ( w — a\\ a 
lim rw + A log = rw H 

ai^a \ w — a J w — a 

with A = a /(a — ai). All logarithmic, rational and 
"mixed" solutions belong to the same class character- 
ized by the property that the first derivative of S + is a 
rational function (or, equivalently, the first derivative of 
S is a mesomorphic function in D(t)). 



T 



In a similar way, one can also consider "mixed" solu- 
tions with both poles and logarithms. In this case the 
function S+ is 



JV K„ 



Ev + E^ 1 ^ 

k—1 rn—1 v / 



(36) 



IV. MULTI-CUT SOLUTIONS WITH 
ANALYTIC CAUCHY DENSITIES OF GENERAL 
TYPE 

A. Motivation 

It is known that all rational solutions (except for those 
of the form f(w,t) = r(t)w + u(t) + A(t)/w, which de- 
scribe a self-similar growth of ellipse) cease to exist in a 
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finite time, because the dynamics gives rise to a cusp-like 
singularity of the interface. These physically meaningless 
singularities just signify that the surface tension effects 
can not be neglected in a vicinity of highly curved parts 
of the interface. Nevertheless, a considerable subclass 
of the purely logarithmic solutions is well defined for all 
positive times and describes a non- sing ular interface dy- 
namics at zero surface tension [l(| I 111. H2I Il3l . [25j . Fur- 



thermore, each logarithmic term in ([31]) has a clear geo- 
metric interpretation of a fjord of oil with parallel walls 
left behind the advancing interface. Its vertex (the stag- 
nation point) is located at Zk — Afclog2, the width is 
7r|Afc|, and the angle between its central line and the real 
axis in the physical plane is argAfc (for more details see 
[T2I . IT3I . I25I]). This interpretation was found to be in an 
excellent agreement with some experiments (see FIG. 2 
in [26| and numerical work (27|). However, more often 
than not, fjords of oil left behind the moving fronts have 
non-parallel walls. Moreover, their walls are not always 
straight, but more often curved, and a non-zero opening 
angle along the fjord was observed [H, [2^] . Such shapes 
of fjords can not be explained by conformal maps with 
a finite number of logarithmic or rational terms. This 
was a significant motivation for us to search for a more 
general class of LG solutions. From mathematical point 
of view, it also looks quite natural to extend the method 
developed above to solutions with singularities of more 
general type than just poles or logarithmic branch points. 

Below we apply the strategy outlined at the end of 
Section [TT] to derive a closed set of equations for the case 
when S+(z) has branch cuts of general type with analytic 
Cauchy densities outside the interface. We start from a 
multi-cut ansatz for S+(z) and then derive an integral 
equation for the conformal map f(w). 



B. A Multi-Cut Ansatz for S+{z) 

The function S-\-(z,t) for z in D{t) is defined as an 
analytic continuation from the air domain, where it is 
given by the integral of Cauchy type ([6]). In the pro- 
cess of analytic continuation one necessarily encounters 
singularities. Typical singularities are branch points and 
poles. In order to keep the function S+(z) single- valued 
one has to introduce a system of branch cuts. The an- 
alytic continuation is achieved by a deformation of the 
integration contour. Moving it towards infinity as far as 
possible, we can write: 



Ei zS+(t)oIt 
res' 



poles 



t(t — z) 



z[5 + (r)] cut dr 
t(t — z) 2ni 

zS+(t) dr 

\t\=R,R-^oo T ( r ~ z ) 2ni 



where [S , + (r)] cu t denotes a discontinuity (a jump) of the 
function [5+(t)] across a branch cut. Let us denote 
branch points by z m , by analogy with logarithmic branch 
points from the previous section. As in that case, and 



for the same reasons, it is convenient to make cuts T m 
from 00 to z m . As before, we assume that the num- 
ber of branch points is finite. Assuming also that S+(z) 
does not have other singularities (in particular, poles) in 
D(t), a single- valued branch of this function can be rep- 
resented as a sum of Cauchy type integrals with Cauchy 
densities V m (j) along the cuts plus a complex constant 
h = S+(0): 



S+(z) 



E 



zV m (r) dr 
r(r — z) 2-ni 



(41) 



The ansate (|4Tj) implies that the jump of S+ on the cut 

r m is 

[S + (z)]=S+{z+)-S + (z-)=V m (z), zeT m (42) 

where z+ (respectively, z_) tends to the point z G T m 
from the left (respectively, from the right) side of the cut 
oriented from the lower limit of integration to the upper 
one. 

There is a big freedom in how to draw the cuts. How- 
ever, if the Cauchy densities can be analytically contin- 
ued from the cuts, then different choices of the cuts with 
the same endpoints are equivalent. Indeed, consider an- 
other system of cuts, T m , connecting the same points. 
Let T m be the Cauchy densities on these cuts. Subtract- 
ing the two integral representations of the same function 
S+(z) one from the other, we see that the "old" and 
"new" cuts combine into closed loops and the Cauchy in- 
tegrals over these loops vanish for all z outside the loops. 
This means that the function V m is just an analytical 
continuation of the function V m , so V m (r) should be re- 
garded as analytic function of the complex variable r 
defined by analytic continuation from the contour L, n . 
Therefore, given analytic functions V m ij) and branch 
points, the choice of cuts T m (if all of them lie in D{t)) is 
irrelevant. For the LG dynamics to be reconstructed from 
the multi-cut ansate (|4"Tj) , all this means that constants 
of motion are the functions V m (z), the branch points z m 
and the constant ti . Note that 00 is a branch point unless 
£ m 7Uz) = 0. 

In this section, we consider only Cauchy densities that 
are analytically extendable without singularities from 
each cut to the whole domain D(t) including infinity. For 
example, one may keep in mind functions regular in D(t) 
including infinity with fixed singularities outside D(t). If 
the Cauchy densities are allowed to have singularities in 
D(t), then the analysis becomes substantially more com- 
plicated. In the next section, we extend the construction 
to the simplest possible singularity of the Vm, a simple 
pole at infinity. This class of solutions includes the im- 
portant case of linear densities. 

It is clear that for Cauchy densities which are constant 
in the physical plane, V m (z) = 2iriA m , the function S+ 
is a linear combination of logarithms. So, the logarithmic 
solutions form a subset of the general multi-cut ones. On 
the other hand, the integral representation (|4"Tj) suggests 
that the multi-cut case can be formally thought of as a 



8 



limit of either rational or (after integrating by parts in 
(|4ip ) multi-logarithmic one, with infinitely many singu- 
larities forming a dense set of points concentrated along 
contours. This remark might be useful for the physical 
interpretation of the multi-cut solutions. 



C. An Integral Equation for the Conformal Map 

The general form of the integral equation obeyed by the 
conformal map f(w) is (fT5|) (or (HOI))- Similarly to the 
rational and logarithmic examples, it can be simplified 
by shrinking the integration contour to the cuts. Namely, 
plugging 



zP m {r) dr 
t(t — z) 2ni 



into the r.h.s. of (115]) . we get 

f(w) — rw + it + 
l_y- I d( [°° f(l/QV m (r) dr 



2ttz ^ /j f , =1 w - C J- Zm t(t - /(1/C)) 2ni 
= rw + uq — 

!_ f°° drV m {j) I /(VC)dC 

t^2 2^ 



(2nif 



C|=i (w-C)(/(l/C)-r) 



Since r lies in the conformal region D(t), there is a unique 
C* inside the unit disk such that /(1/C*) = T - The last 
integral can be calculated by taking the residue at the C* 
(recall that w is outside). After that it is convenient to 
change the integration variable r — > £ connected with r 
by the relation /(1/C) = r. As a result, we obtain the 
following integral equation for the conformal map [37| : 

/(«,) = rw + no - ^ ^ ^ _— . (43) 

Here the points such that z m = /(l/a m ). The 

integration contours 7™ are determined by the change of 
the integration variable. Specifically, consider a contour 
7 m from 00 to l/a m in the mathematical plane, which 
is the pre-image of T m under the map /: T, n = /(7 m ), 
then 7 m is the contour connecting and a m obtained 
as the inversion of the 7 m with respect to the unit cir- 
cle (the transformation w — > l/u>). All the parameters 
in (|43|) except for coefficients of the functions V m (z) de- 
pend on t. The integration contours depend on t as well. 
However, as is explained in the previous subsection, the 
integral does not depend on a particular shape of the 
contours provided they do not intersect each other and 
the boundary of the domain D(t). In what follows we 
do not indicate the integration contours explicitly. Note 
that equation ([4*3")) can be analytically continued inside 
the unit disk provided w does not intersect the cuts from 
to a m . 

In the logarithmic case, V m {z) = 2m A m , the r.h.s. of 
equation (|4"3")) immediately gives the logarithmic ansate 



Y 


k ® 

MS 




Y/w t X 



z=f(w) 



> 




FIG. 1: The correspondence between cuts in the mathemati- 
cal and physical plane. 



for / with moving branch points and constant coeffi- 
cients. In fact it is the integral equation (flB")) that justifies 
the logarithmic ansate. 

Because of importance of the integral equation (|43|) 
we give here an alternative derivation "by hands" which 
is longer but less formal and probably more instructive. 
Using the same arguments as in Section UTTl one can see 
that for each branch point z m of the function S+(z) there 
is a corresponding time-dependent branch point a m of the 
function f(w) determined by the relation z m = f(l/a m ). 
The conformal map can then be written as 



fH 



u 



E 



1 Pm(C) d( 
o 7m C — w 27TI ' 



(44) 



where p m (C) are time-dependent Cauchy densities on the 
cuts 7 m between and o m . 

Let us fix a path T m from 00 to z m lying entirely inside 
the domain D{t) in the physical plane. It is an image of 
a path 7* n connecting the points 00 and 1 / a m outside the 
unit disk: r m = /(7 m ). Let us reflect 7 m with respect 
to the unit circle (i.e., make the transformation w — > 
l/w). The reflected path, j m , connects the points and 
a m inside the unit disk. We take these 7 m to be the 
integration paths in (f44|) . 

Take a point w S 7 TO and consider two points, w+ 
and W- , which are limits to w from respectively left and 
right sides of the branch cut j m (oriented from to a m ). 
Similarly to (|4"2"|) . we have f(w+) — f(w-) = p m (w). The 
corresponding points in the physical plane, z±, tend from 
opposite sides to the point z = f(w*) € T m , where w* = 
1 jw € 7* n is the image of w under the inversion (see Fig. 
[I}. Since conformal maps preserve orientation, we can 
write z± = f(w±). However, the inversion interchanges 
the sides, i.e., w± = l/w T , so z± = f(l/w T ). According 
to equations (fl0|) rewritten as 



S = f(w) 
z = f(l/w) , 



(45) 



the values of the Schwarz function at z± are S(z±) = 
f(wzfi) and the discontinuity of the Schwarz function 
is thus [S(z)} = S(z+) - S(z-) = f(w-) - f(w+) = 
—p m (w). As singularities of S+ and S in D(t) are the 
same, this is exactly the discontinuity of the function 
S + (z). We thus conclude that on j m it holds 



Pm(w) = -V m {z), Z = f(l/w) 



(46) 
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Furthermore, since the r.h.s. is analytically extendable 
to the whole D(t) and the function f(l/w) sets up a 
conformal equivalence between D(t) and the unit disk, 
the p m (w) is analytically extendable from the cut "/ m 
to the whole unit disk. Therefore, (|46|) actually holds 
everywhere in the unit disk. Plugging it into (|44j) . we 
obtain the non-linear integral equation (|43|) for f(w,t). 



D. The Solution Scheme 

To summarize, we have the following formulas which 
express constants of motion through the time-dependent 
conformal map: 



U = 



f(l/w)df(w) 



1 



2tti J\ w \ =1 f(w) 



(47) 



z m = f(l/a m {t),t) . 

They look exactly like the corresponding formulas 
for the logarithmic case. The main difference is that in 
general any explicit representation of the function / is not 
a priori available. It is defined implicitly via the integral 
equation |43|) . 

In principle, these equations supplemented by a for- 
mula for t below give a solution to the problem in an im- 
plicit form, like in rational or logarithmic cases. However, 
the general multi-cut construction is substantially more 
complicated because before applying the scheme outlined 
in section IIIII one has to find a solution to the integral 
equation (|43j) with required analytic properties. Clearly, 
such a solution depends on a mi u and r as parameters: 
f(w) = f(w\r,u , {a m }). Then (|4T|) is a system of equa- 
tions to determine them through the constants of motion. 
The missing equation which includes t can be obtained 
from or ij^]). One of its forms is 



-E 



2ni 



Another form is 



V m (f (I/O) d(f (1/0 -r/C), (48) 



v r ^ 



l.m 



dw pi(Qp m (w) 
2ni (1 - C^) 2 



(49) 



where p m (0 = — ^m(/(l/C))- It 1S an analog of ([55]) 
(compare also with (HU)). Meanwhile, a formula for t\ 
similar to (SHI) also exists: 



P m (f (I/O) dlog((f (I/O), (50) 



which follows from (TlT)) at j = 0. In the case of con- 
stant densities, formulae (|48p and (jST))) immediately give 
expressions (|53"|) and (|32[) for t and t\ from section UlI Bl 



V. MULTI-CUT SOLUTIONS WITH LINEAR 
CAUCHY DENSITIES 

The aim of this section is to elaborate the important 
case of linear Cauchy densities of the Schwarz function 
in the physical plane. This class of solutions is char- 
acterized by the property that the second derivative of 
S+ is a rational function (or, equivalcntly, the second 
derivative of S is a mesomorphic function in D(t)). The 
integral equation for conformal map which is in general 
non-linear, becomes linear in this case. In the context of 
the universal Whitham hierarchy, solutions of this type 
were discussed in [24| (section 7.2) and in [23[ (section 
3.8). This class also includes logarithms and poles on 
the background of the multi-cut functions. However, the 
approach of the previous section is not directly applicable 
because the integrals become divergent. 

These divergences are artificial and can be curbed for 
the price of introducing one more time-dependent pa- 
rameter. The integral equation should be modified. We 
start, in subsection IV Al with a more general situation 
when the Cauchy densities V m (z) are analytic every- 
where in D(t) except for a pole at infinity. In subsec- 
tion IV Bl we specify the results to the most important 
case of purely linear Cauchy densities. In subsection IV CI 
we consider, as an example, some special "finger" pat- 
terns with Zjv rotational symmetry (which are equiva- 
lent to solutions in a wedge with angle 2tv/N). Their 
asymptotic form at large t is given by the known family 
of Zjv-symmetric self-similar "fingers" described by hy- 
per geom etric solutions to the integral equation (|4"3"|) (see 

MM M M M El ) . 



A. Cauchy Densities with a Pole at Infinity 

Let us consider the case when V m (z) are analytic ev- 
erywhere in D(t) except for a simple pole at infinity. One 
immediately sees that neither the multi-cut ansatc (|41[) 
nor the integral equation l[33|) can be directly applied to 
this case because the integrals diverge. Nevertheless, it is 
possible to modify these formulas in such a way that the 
divergences disappear. To this end, consider a modified 
multi-cut ansate: 



S+(z) =h +2t 2 z + 



z 2 V m (T) dr 



oo,r r 



t 2 (t — z) 2iri 



^ (51) 



where t\ and t% are arbitrary complex constants (the first 
and the second harmonic moments) and we adopt the 
same conventions about the cuts T m as before. A single- 
valued branch of this function is fixed by the conditions 
5 + (0) = ti, S' + (0) = 2£ 2 - Let us plug it into the integral 
equation (|19p with d = 1. A simple calculation similar to 
the one done in section HV CI yields the integral equation 



f(w) = rw + u Q + 



in 



E 



(P m (f (l/Q) d( 

w(£ — w) 2iri 



(52) 
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The general scheme of solution remains the same, with 
the only difference that now we have one more time- 
dependent parameter (the coefficient u±). Accordingly, 
we need an extra equation connecting it with integrals of 
motion. The necessary equations, which generalize (|33|) 
and l[32|). can be obtained from fT7|). ([18]) . They are: 



Formulae (1121) read 



UlU Q 



1 C 



/(i/O 



r 



ii = u - 

Y~ E f™ ^(/(VO)d (log(C/(l/C)) - — 
^ E jf p m(f(l/0) d (/(1/C) - £ - «ic) • (53) 



Let us briefly comment on a more general case of 
Cauchy densities with a higher pole at infinity (and ana- 
lytic everywhere else in D(t)). It should be already clear 
how to proceed. If the leading term of V m (z) as z — ► oo 
is z d , then one should extract from S+(z) a polynomial 
Sfcii ktkZ k ~ 1 of degree d and represent the remaining 
part (which is of order 0(z d+1 ) as z — > 0) as integrals 
along the cuts. Plugging this ansate into (fl9|) . one ob- 
tains an integral equation for the conformal map contain- 
ing a Laurels polynomial rw + Ylj=o u jW~i in the r.h.s. 
So, apart from positions of branch points, there are d + 2 
time-dependent parameters r, ito, . . . , 114 which are to be 
connected with integrals of motion by formulae similar 
to (1531). 



B. Linear Cauchy Densities 

The case of linear homogeneous Cauchy densities is 
especially important. Set 



V m {z) = 2mc m z , 



(54) 



where c m are arbitrary complex constants. The explicit 
form of the function S+ is: 



S+{z)=t 1 +2t 2 z + z^c m \og 1 ), (55) 

It is clear that the integral equation (|5"2"|) becomes linear: 



2t 2 = --5><™/ f(l/()d 



1 C 



/(1/C) r 



r U 1 U Q . 
t\ = Uq h 

r 

5>„ r"7(i/od(iog(c/(i/o) 



uqC 



t = r 2 -K| 2 



E g ro y o /(l/o d (/(1/0 - ^ - sic] (57) 

After some simple transformations, they can be brought 
to the form 

Ul 



2t 2 = — 



E 1 



log 













(/(l/C) - r/C) dC 



(58) 



ii = uq — 2u§ti 



E< 



wo log ( - 



rn Zm \ _ 7* 
- MO 



(/(I/O - r/C- So) ^ 
s 

t = r 2 - 2ruii 2 + 
r ^2 

E 



(59) 



rwi log I — 



Y- rU ^2aj: 



P(/(l/C)-r/C-*-«lC)# 

a m JO S> 



(60) 



Let us specify the scheme of solution to the case of lin- 
ear densities. Suppose one is able to find a solution to the 
linear integral equation (|56[) depending on the parame- 
ters a m , u , ui and r: f(w) = f(w\r, uq, Ux, {a m }) (with 
fixed c m ). Then, since f(l/a m ) = z m are constants of 
motion, we get a system of equations for a m , uq, u\ and 
r which becomes closed after adding equations (|58p. (|59[) . 
(|60p . This system determines the parameters as implicit 
functions of t. 

We conclude the subsection with a remark that the 
solutions of the linear integral equation (|56|) for the con- 
formal map f(w,t), which correspond to linear Cauchy 
density of the Schwarz function defined by (|55p . deserve 
a special name in view of their importance. So, in what 
follows we will refer to them as to hyper-logarithmic so- 
lutions, since they contain logarithmic solutions (|31[) and 
the Gauss %per-geometric function (see (|55|) below) as 
particular cases. 



Example: Zjv-Symmetric Solutions with N 
Radial Cuts 



In this subsection we consider perhaps the simplest 
non-trivial application of the linear integral equation de- 
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rived above: growing patterns with Zjv rotational sym- 
metry whose Schwarz function has exactly N radial 
branch cuts in D(t) with linear Cauchy densities. 

Let uj = e 27Tl / N be the primitive root of unity of degree 
N. The Zjy-symmetry implies S(wz) = uj^ 1 S(z). More- 
over, this relation holds for both + and — parts of the 
Schwarz function separately: S±(ujz) = lj^ 1 S±{z). This 
prompts us to choose the constants of motion in the form 



C„, = CUJ 



-2m+2 



= \LU % 



m = 1,2,..., TV, (61) 



with real positive constants c and A, so the function 5+ 
is 



JV-l 



5 + ( Z ) = C z^ W - 2m log 1 



m— 



ZUJ 



-m+1 



A 



cz 



N-l 



(i-4)a w 



_~2 N ,{*/< 
Z N 



where 2-F1 is the Gauss hypergeometric function. 

For the conformal map the Zjv-symmetry means 
f(uj k w) = ui k f{w) 1 so the expansion of the function 
f(w)/w goes in powers of w~~ N : 



f(w)/(rw) = l + J2f< 



kW 



-kN 



k>\ 



where fk = u^n /f- The reality of c and A in (|62[) implies 
that the coefficients fk are real, i.e., f(w) = f{w). In 
what follows we assume that N > 3, so the coefficients uq 
and Mi in the expansion of f(w) ([2]) are always identically 
zero. 

We use the integral equation ([56]) . Substituting our 
data, we get 



JV-l 

f{w) =rw + c u> 

m=0 



2 m 



w(( — w) 



(63) 



where a is a point between and 1 such that /(1/a) = A. 
Choosing the integration paths to be straight lines from 
to aui m and using the Z^-symmetry, we arrive at the 
equation 



/H 



cN 



1/11 



wf(x) dx 
1 - w N x N 



(64) 



where c is a constant and r, a depend on time. Let 
f(w) = f(w\r,a) be a solution to this integral equation 
with parameters r, a, then three parameters r, a and t 
are connected by two equations which determine r, a as 
implicit functions of t. The first of these equations is 
obtained from the fact that A = f(l/a(t),t) is a constant 
of motion by substituting w = 1/a into (1641) : 



= Xa-cN 



f(x\r, a) dx 



1/a 



l-x N a 



N„-N 



(65) 



the second is one or another version of the formula for t 
in terms of f(w). 



UK 



np/2 



air 



FIG. 2: Laplacian Growth in wedge geometry. Air is pushed 
into oil trough the corner of the wedge. f3 is the fjord angle. 



We have obtained a one-parametric family of exact so- 
lutions to the LG describing growth of N symmetric air 
fingers or, what is mathematically the same, growth of 
an air finger in the wedge with interior angle 2tt/N (Fig. 
[2]). It is not clear at the moment whether the conformal 
map f(w) is available in a closed analytic form. However, 
the solution can be analyzed numerically (see appendix 
©. 

What is the meaning of the parameter c? To answer 
this question, consider the limit t — > 00, wherein, as one 
can easily show, a — > 1. Then equation ([M]) becomes 
identical to the linear integral equation describing self- 
similar growth of Zjv-symmetric fingers (or fingers in the 
wedge with interior angle 2t:/N): 



/M 



N sin 7r/3 



wf{x) dx 
1 - w N x N 



(66) 



with the "boundary condition" f(l) = 0. As shown in 
our earlier paper [16| . (3 is the interior angle of the oil 
fjord between two neighboring air fingers. Since c in 
is a constant of motion, we conclude that setting 



sin it (3 



(67) 



we can interpret (3 as the asymptotic fjord angle. 

Meanwhile, an analytic solution of (f6"6"|) in terms of the 
Gauss hypergeometric function is available: 



f(w) = rw{l - w~ N f 2 Fi 



/3, 0-K 
1 - K 



-N 



(68) 



where 



For more details, see [1| 
studied in detail. 



2 

;= AT 

where the self-similar case was 
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The integral equation (|64|) is equivalent to an infinite 
system of linear equations for the coefficients fk- To rep- 
resent it in this form, let us introduce a function F via 

oo 

f(w)=rwF(w- N ), F(w)=J2fkW k , /«, = 1 , 

then the integral equation (|64[) becomes 

r./ \ 1 f a x- K F(x)dx 

F(w) = 1 + c / , (69) 

where a = a N . Substituting the Taylor expansion of F 
and comparing the coefficients, we get 

/i = -"Ej+F^^' (70) 

(However, because the r.h.s. contains /o = 1, this is 
not an eigenvalue equation but rather an inhomogeneous 
system of linear equations.) It is also useful to note the 
formula for t, 

t = r 2 + iiVcA 2 + cr 2 V !± , (71) 

2 — ' n — k v ' 

n>0 

which follows from (f^Oj) . 

VI. DISCUSSION AND CONCLUSION 

In conclusion, let us briefly comment on the following 
three major points of this work: 

• The inverse potential problem integral equation, 
03} (or (20])), 

• The equations for multi-cut solutions, (|41 [) and 

• The equations for hyper-logarithmic solutions, (|55p 
and ijgB]), 

The inverse potential problem equation (|15p : The formu- 
lation of the LG problem as a linear growth of a domain 
area with conserved harmonic moments requires to solve 
the inverse potential problem, namely to find the confor- 
mal map f(w) from the function S+(z) = Y^kLi ktkz k ~ 1 , 
whose Taylor coefficients are conserved harmonic mo- 
ments of the growing domain and from the area/time 
t. While it is straightforward to obtain S+(z) and t from 
/(u>) (the direct potential problem), it is a well-known 
challenge to do it the other way around (the inverse po- 
tential problem) The reformulation of the inverse 
potential problem as the nonlinear integral equation (|15|) 
(or pp)) ) for the function f(w), assuming that S+ (z) is 
known, appears to be extremely useful. In particular, it 
allows us to obtain some important classes of solutions to 



the LG problem with prescribed integrals of motion (har- 
monic moments), including those with arbitrary branch 
cuts as well as those with pole singularities. This inte- 
gral equation is also expected to be of significant help in 
a future work on LG and related problems of interface 
dynamics. 

The equation for multi-cut solutions, (|43[) : A mathe- 
matical motivation for this work was a strong feeling that 
rational and logarithmic conformal maps do not exhaust 
the list of exact solutions of LG at zero surface tension. 
A physical motivation came from a long standing need to 
describe most general moving interfaces in Hele-Shaw ex- 
periments with negligible surface tension. As was already 
said in section IIV Al while finite linear combinations of 
logarithms can describe interface dynamics without finite 
time singularities, they fail to describe an interface with 
non-parallel fjords walls. This motivated us to look for 
a more general family of solutions to the LG equation 
©. We expect that the multi-cut solutions presented in 
this paper do explain the experiments [28|, l29l . [36[ where 
formation and development of oil fjords with non-parallel 
and/or non-straight walls was observed. 

Two special kinds of branch points and cut singular- 
ities were already considered in earlier works on Lapla- 
cian growth. Fractional time-independent branch points 
appeared in studies of self-similar singular interfaces 
[3, EE EH , while logarithmic branch cuts with constant 
Cauchy densities appeared in [l(], EH E3, El, l25j . We 
have shown that all of them are just various special cases 
and limits of the presented general construction. In the 
light of our approach, it also becomes clear why the pole 
and logarithmic solutions are special: the Cauchy densi- 
ties are particularly simple in these cases, so the integral 
equation can be easily solved. It also shows a road to 
obtain new families of exact solutions. 

The equation for hyper-logarithmic solutions, (I5GD : 
The mathematical significance of the hyper-logarithmic 
solutions is that they correspond to the linear integral 
equation ([56]) . which ought to be much more accessible 
to analytic treatment than the general nonlinear equa- 
tion (|43[) for multi-cut solutions. The physical impor- 
tance lies in the belief that the hyper-logarithmic solu- 
tions describe fjords with a constant opening angle, in 
agreement with viscous fingering experiments [29| . The 
interface dynamics simulated in a wedge, shown on Fig- 
ure 02 seems to confirm this belief, if to consider fjords 
central lines as walls of virtual wedges in accordance with 
[28| . A thorough geometric analysis of the corresponding 
interface dynamics based on the integral equation (|56p 
will be published elsewhere. 
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APPENDIX A: NUMERICAL SOLUTION OF 
EQUATION ([64]) 

Equation (|64|) with c given by (|67p can be rewritten in 
the following way: 

N sin 7r/3 / /(C,*)log(l-a(t)C) dC 
f(w,t)=r(t)w —*>f ^ w n C n Sri- 



f(l/a(t),t) = X, 



(Al) 



We are looking for a solution consistent with the Zjy- 
symmetry: 



f(w,t)=r(t)w + J2f k (t)i 



,-Nk+l 



(A2) 



fc=i 



Equation (IA1I) then takes the form 

sin 71-/3 a n - K (t) 



/»(*) = -r(t) 

7Vsin7r/3 ^ a 
>. Jk[t)- 



7T n — K 

n-\-k — K 



it) 



(A3) 



fc=i 



where k = 2/N, a = a N . Using the relation A = 
f(l/a(t),t) and re-denoting f k (a) = f k (t(a)), r(a) = 



r(t(a)), we get 



■(a) =a K/2 A-^/ fc (a)a fc , 



(A4) 



fe=i 



and 



E' , . fca" +fe 7r(n — k) 



k=l 



71 + k — K Sin 7T/9 



This equation is easily solved numerically for any given 
a . The conformal radius r and the time are then found 
from (|A4|) and (J22J) (or f?T|)) respectively. The solution 
is shown in figure [3J 




0.5 



FIG. 3: (color online) The results of the numerical solution 
of the equation ([64]) for N = 6, /3 = 0A/N. The left panel 
shows the dependence of the area t on the conformal radius 
r, the right panel shows the growing finger at different times. 
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